Two Bounds on the Maximum Phonon-Mediated Superconducting Transition 
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Two simple bounds on the T c of conventional, phonon-mediated superconductors are derived 
within the framework of Eliashberg theory in the strong coupling regime. The first bound is 
set by the total electron-phonon coupling available within a material given the hypothetical abil- 
ity to arbitrarily dope the material. This bound is studied by deriving a generalization of the 
McMillan- Hopfield parameter, rj(E), which measures the strength of electron-phonon coupling in- 
cluding anisotropy effects and rigid-band doping of the Fermi level to E. The second bound is set 
by the softening of phonons to instability due to strong electron-phonon coupling with electrons at 
the Fermi level. We apply these bounds to some covalent superconductors including MgB2, where 
T c reaches the first bound, and boron-doped diamond, which is far from its bounds. 

PACS numbers: 74.20.-z, 74.62. Dh,74.62.Bf,74.25.Kc 



I. INTRODUCTION 

With the discovery of increasingly high transition tem- 
peratures for conventional BCS superconductors, it is of 
interest to reassess and rederive bounds on the maxi- 
mum possible T c . The goal in this area of research is 
to provide useful theoretical guidelines for the continued 
search for higher T c materials. In the work of Cohen and 
Andersoni, three different arguments give three different 
bounds on T c . The first of these arguments is simply 
based on the form of the McMillan T c formula 2 , which 
gives a maximum T c proportional to a material's average 
phonon energy oj p h for all values of the electron-phonon 
coupling constant A. However, the McMillan formula is 
only valid in the A < 1 regime and the correct Eliashberg 
T c formula for large A is^ 



(1) 



which is unbounded in A and potentially much larger 
than u)ph- The second argument of Cohen and Ander- 
son postulates that stability requires the static dielectric 
function of a material to be positive, which leads to an 
upper bound on T c of approximately 10 K. This bound 
and stability requirement are subsequently dismissed in 
the paper because of the neglect of local fields and umk- 
lapp processes, and further work* has clarified that sta- 
bility only requires a positive static dielectric function at 
long wavelengths. This second argument is refined and 
it is concluded that strong electron-phonon coupling is 
limited by an instability leading to insulating, covalently 
bonded materials. Strong covalent bonds and supercon- 
ductivity are not mutually exclusive, most notably in the 
case of MgB^, but engineering their coexistence is a dif- 
ficult materials problem. The final argument examines 
the structural instability resulting from phonon softening 
due to the same electron-phonon coupling that causes su- 
perconductivity. Within the framework of a simple semi- 
conductor model, it is shown that a phonon can go un- 
stable before T c reaches a maximum, further reducing 



the bound from the naive T c ~ u p h . Expressions relating 
phonon softening to A have since been generalized£iii2i2ii£ 
and take the form 



(2) 



where Q p h is a bare phonon frequency. 

These previous studies are revisited here in a more 
general framework to attempt to obtain more rigorous 
upper bounds on T c . First, the strong coupling limit 
of Eliashberg theory is clarified through the use of up- 
per and lower bounds on T c . An approximation for T c 
and electron-phonon coupling are directly related to the 
properties of the constituent atoms and their local envi- 
ronments within a material. We examine this measure 
of electron-phonon coupling as a function of electron en- 
ergy to give insight into how a material might achieve a 
higher T c based on modifications of its electronic struc- 
ture. We then derive a structural instability bound on T c 
through the calculation of the phonon softening caused 
by the same phonon-mediated pairing interaction that 
forms Cooper pairs near the Fermi surface. These meth- 
ods are applied to probing the limits of superconductivity 
in three boron-carbon covalent metals: MgB2, Lii-^BC, 
and boron-doped diamond. 

The primary result of this paper is that there are two 
fundamental bounds on the transition temperature of a 
phonon-mediated superconductor, one set by lattice sta- 
bility and the other by the total strength of available 
electron-phonon coupling. The lattice stability T c bound 
can be approximately derived just from equations Q and 

. Through substitution, the strong coupled T c formula 
can be written without A as 



9. 
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(3) 



which is maximized as the electron-phonon coupling 
grows, A — > co, softening the phonon frequency towards 



2 



instability, w p /, — > 0. This suggests that Eliashberg the- 
ory does not have an unbounded T c as Eq. (JIJ might 
lead us to believe, but rather it is proportional to a bare 
phonon frequency £l p h that can be larger than the soft- 
ened phonon frequencies observed in a superconducting 
material. Another form that Eq. JQ) is sometimes writ- 
ten as is 



a 

M' 



(4) 



where M is an averaged atomic mass and r] is the 
McMillan-Hopfield parameteA The rj has units of a 
spring constant and is related to the strength of the elec- 
tronic response of electrons at the Fermi surface to per- 
turbations of the atoms in the crystal. For a given M, 
a T c bound is set by the maximum rj possible consid- 
ering all energies at the Fermi level accessible through 
chemical doping or other physical modifications of the 
material. A large rj is required to approach the T c bound 
set by phonon stability in materials with strong covalent 
bonds and a large £l p h- Approaching the phonon stability 
bound in systems containing strong covalent bonds will 
require careful material design. 



II. ELIASHBERG T c BOUNDS 

The T c calculated from the isotropic Eliashberg equa- 
tions has been the subject of both upper and lower bound 
studie d 11 " . In the isotropic case, the T c of a material 
is uniquely determined by the Eliashberg spectral func- 
tion a 2 F(ui) that characterizes the electron-phonon cou- 
pling and the Coulomb parameter /i* that quantifies the 
strength of the screened electron-electron Coulomb re- 
pulsion. Approximate formulas for T c usually represent 
the average electron-phonon coupling by 



A = 2 / a 2 F{uj)u)- x duj, 
Jo 



(5) 



and an average phonon frequency defined by a weighting 
function f(u>), 

2 f°° 

(f( u )) = -J f(aj)a 2 F(uj)uj~ 1 duj. (6) 

At large temperatures, T only appears in the Eliashberg 
equations as X(lo 2 )/T 2 , which is where the asymptotic T c 
formula of Eq. (JJJ originates^. If only X(lo 2 ) is specified, 
T c cannot be precisely determined, but it can be bounded 
from above and below, 



< T c < 0.183y/X(uj 2 ) 



(7) 



Both bounds can be attained by an Einstein model 
phonon spectrum, 



the lower bound corresponding to loe — * oo and the upper 
bound corresponding to lo e ^ C&ii. The spectra associ- 
ated with the two bounds are both unphysical patholo- 
gies, corresponding to the cases of maximum phonon soft- 
ening and infinite phonon hardening. 

We attempt to tighten the bounds by independently 
constraining A and (lo 2 ) as well as specifying fi* . It is 
already known^ that the Einstein spectrum maximizes 
T c for a given A and (lo 2 ). However, the lower bound 
on T c is still zero even with these additional constraints. 
This is due to a pathological spectral function with two 
peaks, one approaching zero and the other diverging to 
infinity. In order to obtain a nontrivial lower bound, 
additional constraints need to be placed on the spectral 
function, and we choose the added restriction of a max- 
imum phonon frequency Lo max that replaces the upper 
limit of integration in Eq. JSJ and ©. The spectral 
function that minimizes T c at fixed A, (lo 2 ), and u> max is 
numerically determined^ to be 



a 2 F{to) r 



Xcv 
2U 2 

max 



(( 



(lo 2 ))S(lo)+ 

(uJ 2 )S(lO - Umax))- (9) 



The phonon at lo — contributes to A but not to T c , 
making this effectively an Einstein spectrum when calcu- 
lating T c . The qualitative result here is that T c is maxi- 
mized when the spectral function has no width and it is 
minimized when all the spectral weight is pushed to the 
edges of the allowed frequency interval. The stronger, 
more complicated T c bounds take the form 



/ (X(L0 2 )/L0 2 max ^*) Lo max <T c <f (A, //) vV). (10) 

The function /(A,/x*) smoothly interpolates between the 
exp(— 1/A) weak coupling limit and the y/~X strong cou- 
pling limit and can be empirically fit to 



f(X,n*) = 0.69 exp 




(11) 



a 2 F(oo) e — Xloe5(lo — loe) /2, 



(8) 



for A > fJ,*, otherwise / = 0. This functional form is 
chosen to balance simplicity and accuracy and to match 
exact analytic results— i for < A < 0.4 and A> 1. The 
coefficients were chosen based on a least squares fit on 
the interval 0.4 < A < 20 and < ji* < 0.2. The average 
root-mean-squared deviation of / is 0.01 and the max 
deviation is 0.03. 

We plot the upper and lower bou nds o f T c for a rea- 
sonable but arbitrary uo max = 1.5^/ (lo 2 ) as well as the 
McMillan formula for an Einstein spectrum in figure ^ 
The McMillan formula deviates perceptibly from the up- 
per bound around A = 0.5, by A = 1 the deviation is 
roughly 20%, and after A = 2 it passes below the lower 
bound chosen for the figure. The relative width of the T c 
bound, (T™ ax -T™ ln )/T™ ax , is largest at weak coupling. 
It has been observed that T c is more accurately propor- 
tional to a logarithmically averaged phonon frequency, 
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FIG. 1: Comparison of T c upper and lower bounds (assuming 
Umax = l.5y/ (lo 2 )) with the McMillan formula for an Einstein 
spectrum of frequency Lu p h with a varying electron-phonon 
coupling strength A and fi* = 0. 



exp((ln(w))), rather than Lo max or yj (lo 2 ) in the weak- 
coupling regime^. As a result of that analysis, the most 
widely used formula for calculating T c is the McMillan 
formula modified to include a log-averaged phonon fre- 
quency. In the strong coupling regime that is of more 
interest here, A > 1, both upper and lower bounds on T c 
are proportional to */ X(lo 2 ) and converge as A grows. 

To complete our discussion of T c bounds, we need to 
address the effects of anisotropy. The anisotropy of the 
Eliashberg spectral function causes an increase in T c from 
that calculated in the isotropic approximation 14 . The 
reason for this enhancement is an increase in A due to 
anisotropy, most notably essential for describing the su- 
perconductivity in MgBa 1 ^* 1 ^, raising A from 0.77 to 1.01. 
The simplest model in which to consider this effect is a 
multiband superconductor, which introduces a band in- 
dex to the Eliashberg equations!!. There is a separate 
gap value for each band i, A{(cj), and A and X(lo 2 ) are 
separated into phonon processes that couple band i to 
band j, 



and 



Ajj = 2 



[A(w 2 )]j 



a 2 Fij(u>)u) 1 cLo 



a 2 Fij(u>)ucLj. 



(12) 



As in the isotropic case, for high temperatures the Eliash- 
berg equations only contain T in terms of the form 
[X(u 2 )]ij /T 2 . Correspondingly, the simple bounds on T c 
are the same as Eq. (J7J with X(cu 2 ) replaced by the largest 



eigenvalue of [X(u> 2 }]ij. Both A^ and [A(w 2 )] y 
lated back to the isotropic approximation by 



□ 



E 



Nj(E F ) 
N(E F ) 



can be re- 



(14) 



where □ is A or X(lo 2 ), Ni(E F ) is the density of states in 
the i th band at the Fermi level, and N(E F ) = J2i Ni(E F ) 
is the total density of states at the Fermi level. The 
largest eigenvalue of Dy can be written variationally as 



^XjXjNjjEF)^ 



(15) 



The isotropic value corresponds to Xi = 1, less than or 
equal to the maximum value. Starting from the isotropic 
lower bound spectral function in Eq. @ , any anisotropic 
perturbations at fixed A, (lo 2 ), and Lo max can only in- 
crease the largest eigenvalues of Ay and [X(co 2 )]ij and 
thus only increase T c . The lower bound of Eq. 110|) 
cannot be lowered by anisotropy and still applies in the 
anisotropic case. For simple calculations of T c upper 
bounds, we define effective electron-phonon coupling ex- 
pressions that don't require detailed information about 
the anisotropy of a system, 



{LO 



A = y Xu 

i 

^A-*]T[A<o 



(16) 



Both A and X(uj 2 ) will always be larger than the largest 
eigenvalues of Ay and [X(to 2 )]ij respectively because they 
are each the trace of a positive definite matrix. In the 
isotropic case, A equals A and (lo 2 ) equals (lo 2 ). Starting 
from the isotropic upper bound spectral function in Eq. 
©, any anisotropic perturbations at fixed A and (oj 2 ) can 
only redistribute the magnitude of the largest eigenvalues 
of Xij and [A(w 2 )]y to smaller eigenvalues and thus only 
decrease T c . Another important effect of anisotropy is in 
reducing the effects of //* on the effective fi* being 
determined by the distribution of the gap between bands, 



»eff 



(17) 



which in cases of extreme anisotropy can completely nul- 
lify the effect of fx* . The end result is that the upper 
bound of Eq. (|10fl still applies to the anisotropic case 
provided that A and (w 2 ) replace the isotropic values and 
fi* is set to zero. 

The anisotropic case now provides additional difficulty 



(13) in the strong coupling limit because X(lo 2 ) and A( 



don't necessarily converge to each other and there is 
no unique simple measure of T c that is independent of 
anisotropic details. Rather than resort to a more compli- 
cated expression, we examine the accuracy of A and A in 
a common anisotropic scenario. Most conventional BCS 
superconductors are reasonably isotropic and one of the 
most anisotropic BCS superconductors known, MgB 2 , 
fits well into a two band modeli^i&. In this model, there 
is a strongly coupled cr-band Fermi sheet and a more 
weakly coupled 7r-band sheet. We take a more extreme 



4 



version of this two band case, where the second band has 
no electron-phonon coupling at the Fermi surface. In this 
case, the correct anisotropic result is [A(w 2 )]n, which is 
also equal to the upper bound X(oj 2 }. The isotropic ap- 
proximation X(lu 2 ) equals Ni(Ef)/N(Ef)[X(uj 2 }]ii and 
can substantially underestimate the correct result de- 
pending on how much of the density of states at the 
Fermi level lies in the decoupled band. In general, A will 
be correct for single-gap superconductors and the error 
for multi-gap superconductors will be proportional to the 
total magnitude of all gaps excluding the largest. For this 
reason, we consider X(ui 2 ) to be a better approximation 
than A(ct> 2 ) for calculations of anisotropic superconduc- 
tivity. 

Assuming that we are able to enter the strong coupling 
regime, we are guaranteed to increase T c uniformly if we 
focus solely on increasing A(cD 2 ). This approach works 
especially well if we consider classes of materials with 
similar Q, p h because Eq. @ shows that increasing X(ui 2 ) 
also increases A due to phonon softening. It can also 
be detrimental to independently consider A and (uj 2 ) in 
classes of materials such as is portrayed in Eq. 10. In- 
creasing the coupling to the very low frequency modes 
of a system most efficiently increases A but doesn't raise 
T c . In such an extreme case, it is more useful to have 
a metric for superconductivity that is insensitive to low 
frequency coupling, which is true of X(uj 2 ). We make no 
attempts to model /i* as it is typically between 0.1 and 
0.15 in metalsiS and variations of this magnitude can 
only change T c by a few percent in the strong coupling 
regime. We are left with a single quantity to optimize, 
A(w 2 ), that serves as a simple metric for enhancing T c . 



III. LOCAL APPROXIMATIONS 

As was first pointed out by McMillan^, the quantity 
X(uj 2 ) is independent of the details of the phonon spec- 
trum. It can be written as a sum of atomic contributions 



x(uj 2 ) = — , 

v 1 ^ M, ' 



(18) 



where i is summed over all atoms in the unit cell, Mi 
are the atomic masses, and rji are the McMillan-Hopficld 
parameters of the atoms. The quantity A(w 2 ) is also inde- 
pendent of the phonon spectrum and can also be written 
as a sum of effective McMillan-Hopfield parameters, 



v ' ^ Mi 



(19) 



From the definition of X(uj 2 ) in Eq. (|To|) and McMillan's 
derivation of r\ it follows that 



nk 



where dVscF /dRi is the change in the self-consistent 
potential due to a change in the ith atomic position, and 
|n,k) and E n k are the Bloch wave function and band 
energy for band n and wavevector k. Just as with A, rji 
is less than or equal to rji with equality occurring in the 
isotropic case. We can generalize the electron-phonon 
coupling parameter rji to electrons away from the Fermi 
surface, at an energy E, simply by substituting Ep — > E, 



i /k 



OVscf 

(n ' k| ^RT |n ' k) 



5(E - Erik). (21) 



5{E F - £ nk ), (20) 



This expression has the form of a density of states 
weighted by the sensitivity of the electronic states to 
atomic displacements. Whereas A = and rji(Ep) = 
for all non-metallic materials, rji(E) will have non-zero 
values and potentially interesting features for all materi- 
als. By looking at this coupling parameter as a function 
of energy, we can immediately examine rigid-band dop- 
ing scenarios and get a better sense of the magnitude of 
hypothetical electron-phonon couplings. 

The effects of anisotropy have been considered in de- 
riving rji, but so far neglected is phonon anharmonicity, 
which can be of equal importance in materials such as 
MgBi 15 . Anharmonicity is important for "soft" phonon 
modes where the average displacement of atoms is large 
enough that it samples the non-harmonic terms of the po- 
tential energy surface. Eliashberg theory can be formally 
reformulated to account for anharmonicity^, where the 
largest effect is typically an increase in phonon frequen- 
cies and other corrections contain multi-phonon processes 
and Debye- Waller factors. The most common approach 
taken is the inclusion of anharmonic phonon frequencies 
wherever phonon frequencies appear in Eliashberg the- 
ory. At this level of correction, rji and thus the quan- 
tity A (w 2 ) are unchanged as they are independent of 
the phonon spectrum. Anharmonic corrections can of 
course shift the average phonon frequency (uj 2 ), but A 
will correspondingly change such that their product re- 
mains constant. Taking into account more anharmonic 
effects, rji can also change because Vscf in Eq. i|2U|) 
is dependent on atomic positions and will be affected by 
large atomic displacements. The anharmonic phonon fre- 
quency correction is significant because as a phonon soft- 
ens the harmonic change in energy with phonon displace- 
ment, cPE/d^ 2 , gets small and the phonon frequency 
becomes dominated by anharmonic corrections propor- 
tional to (d 4 £/d£ 4 )Ai? 2 , where AR is the magnitude of 
the average phonon mode displacement. In contrast, the 
anharmonic rji correction typically will not be important 
because the harmonic value does not get small as phonons 
soften and as is shown in section HVI one mechanism of 
phonon softening is actually related to large rji values. 
The only anharmonic effects we will consider are correc- 
tions to (uj 2 ) and not rji. 

To calculate rji (E) using a modern electronic structure 
code, we must find an approximation to dVscF / '<9R; in 



Eq. Q21[l. The potential due to the displacement of a 
single atom in a crystal breaks translational symmetry 
and is a difficult quantity to compute using codes that 
rely on periodicity. A reasonable approximation to this 
single atom displacement is the displacement of a single 
atom per unit cell within a sufficiently large supercell of 
the material of interest. This displacement is a zone- 
center perturbation and the matrix elements in Eq. i|21|) 
are simply obtained from first order perturbation theory 
on the band energies, 



dE n]s 



dR, 



8{E - E nk ) 



(22) 



This calculation is no more difficult than decomposing 
the force on each atom into its individual band and 
fc-point contributions and then calculating a density of 
states. The rjt obtained from this procedure does not re- 
quire the rigid-mufRn-tin approximation (RMTA) often 
made in other more simplified calculations of electron- 
phonon couplingSiS. By utilizing a supercell approx- 
imation, we assume spatial locality of electron-phonon 
coupling as in the RMTA, but here this assumption can 
be systematically verified and improved by increasing 
the size of the supercell. The calculations in this pa- 
per are performed using the PWSCF electronic structure 
code2£ with the provided default norm-conserving pseu- 
dopotentials and converged planewave cutoffs. Unless 
otherwise noted, calculations are performed with exper- 
imental lattice constants. Density of states and rji(E) 
calculations are performed using gaussian smearing and 
dE n k/dRi values based on multiple calculations with dis- 
placed atomic coordinates and a finite difference approx- 
imation. The gaussian smearing width and k-point grid 
utilized in each calculation are noted with the results. 
The rji(E) values depend on the choice of unit cell, and 
for easier comparison we use a unit cell-independent rj(E) 
value that is summed over the atoms in the unit cell. The 
summed rj(E) value can be used in Eq. (|19|l with an ap- 
propriate average atomic mass. As an example, we calcu- 
late rj(E) for diamond in figure [21 Within the precision 
of the calculation, there is no difference between rj(E) 
calculated in one unit cell and in a 2 x 2 x 2 supercell. 
This demonstrates the locality of the electron-phonon 
coupling that results from a strong covalent bond. Be- 
cause a large supercell isn't required for accurate rj(E)'s, 
this is a computationally efficient methodology for calcu- 
lating electron-phonon couplings in systems of this type. 

Based on this analysis, we can now write T c in the 
strongly coupled limit as 



(23) 



where Ff s is the average root-mean-square force on 
atom i from electrons on the Fermi surface. This re- 
lation reflects the common wisdom based on the BCS 
model that a large density of states at Ep and light ele- 
ments are what lead to a high T c . It also helps to give a 




FIG. 2: Density of states (DOS) (dashed line) and rj(E) (solid 
line) for diamond. This calculation is performed on one unit 
cell of diamond with a 24 x 24 x 24 k-grid and 0.03 Rydberg 




FIG. 3: Density of states (dashed line) and rj(E) (solid line) 
for silicon in the diamond structure. This calculation is per- 
formed on one unit cell with a 24 x 24 x 24 k-grid and 0.015 
Rydberg gaussian smearing. 



simple qualitative picture of electron-phonon coupling by 
allowing us to focus on a single electronic state and atom 
at a time. If the electron occupies a single orbital of this 
atom, its charge distribution will be symmetric about the 
atomic center and it will contribute no net force to the 
atom. This is related to the I — > Z'±l angular momentum 
sum rule in the electron-phonon coupling calculations of 
Hopfield^. Nonzero coupling arises from deviations from 
spherical symmetry around an atom and the length scale 
for this deviation is the inter-atomic spacing. In terms of 
a tight-binding model with inter-atom hopping param- 
eter t, we can use a simple chain rule argument to re- 
late rj to t, dE/dR = (dE/dt)(dt/dR), whereby a atomic 
displacement perturbs the neighboring hopping parame- 
ters and they in turn perturb the band energies. Strong 
electron-phonon coupling is related to a large sensitivity 
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FIG. 4: Density of states (dashed line) and rj(E) (solid line) 
for cubic boron nitride. This calculation is performed on one 
unit cell with a 24 x 24 x 24 k-grid and 0.03 Rydberg gaus- 
sian smearing. In a single unit cell approximation, rji(E) is 
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FIG. 5: Density of states (dashed line) and rj(E) (solid line) 
for graphite. This calculation is performed on one unit cell 
of AA-stacked graphite with a 24 x 24 x 12 k-grid and 0.03 
Rydberg gaussian smearing. Only in-plane motions of the 
carbon atoms are considered; the out of plane motions do not 
contribute significantly to rj(E). 



in t with respect to atomic displacements, dt/dR, and 
not directly to a large t. 

To study the effects of covalent bonding on electron- 
phonon coupling, we compare rj(E) values of several re- 
lated materials. As an example of changing bond length 
and strength, we compare the rj(E) of diamond in hgure 
[5] with silicon in a diamond structure in figure [3] Carbon 
has rj{E) values between two and ten times larger than 
silicon presumably due to a shorter, stronger covalent 
bond and less core screening. The density of states of 
silicon is twice that of carbon and within a tight-binding 
picture this is associated with a factor of two reduction 
in t. Although we generally associate larger A values to 
larger density of states values at the Fermi surface, this 




FIG. 6: Density of states and rj(E) for simple cubic tight 
binding in D=l,2,3 dimensions. 



reduction in t is concurrent with a reduction in dt/dR 
that leads to an overall decrease in rj(E). We can next 
examine the effect of bond ionicity on rj by comparing 
diamond with boron nitride in a zinc-blende structure 
in figure 01 Although their bond lengths are very sim- 
ilar, boron nitride has a more ionic, less covalent bond 
than diamond and the dt/dR and rj values are corre- 
spondingly reduced. All these materials share the dia- 
mond (zinc-blende) structure and contain similar peaks 
in rj(E) above and below the Fermi level, due to strong 
coupling associated with the anti-bonding and bonding 
states respectively. 

Besides changes in bond strength and ionicity, we can 
also change the character of a covalent bond by chang- 
ing the structure of a system. Carbon in graphite has 
sp 2 bonding instead of the sp 3 bonding of the diamond 
structure. In this structure there are only three bonds per 
atom rather than four in diamond, but the bond lengths 
are shorter and the bonds are stronger. In figure [S] we 
plot the density of states and rj(E) of graphite. The rj(E) 
for valence states of diamond and graphite are more sim- 
ilar in shape and magnitude than their density of states. 
The biggest difference in valence elcctron-phonon cou- 
pling between these systems is the peak in coupling of 
diamond being replaced by a smaller, wider plateau in 
graphite. The features of rj(E) do not directly correspond 
to density of states features in the diamond-structure ex- 
amples, and the deviation is even more pronounced in 
graphite. In particular, the maximum in the valence den- 
sity of states of graphite associated with a van Hove sin- 
gularity of graphene is not associated with a maximum 
in electron-phonon coupling. 

To understand the fundamental difference between 
N(E) and rj(E), we calculate both for a simple cubic 
nearest-neighbor tight binding model with no on-site en- 
ergy and a hopping energy t in one, two, and three di- 
mensions. Electron-phonon coupling is modeled by as- 
suming that atomic displacements at a site correspond 
to changes in hoppings from that site parameterized by a 
constant dt/dR dotted with the unit normal vector in the 
direction of the hopping. The N(E) and rj(E) of these 
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FIG. 7: Density of states in the primitive cell (solid line) 
and rj(E) of fee aluminum for a 2 x 2 x 2 supercell with an 
18 x 18 x 18 k-grid (dotted line) and a 3 x 3 x 3 supercell with 
a 12 x 12 x 12 k-grid (dashed line). The experimental value 
of tj at E F is 2 eV/A 2 . 



three systems are plotted in figure El showing that van 
Hove singularities characteristic of low dimensionality do 
not cause divergences of rj(E) . The shape of rj(E) is rel- 
atively insensitive to dimensionality and is proportional 
in all cases to a generic tight binding result, 



dt 



dR 



(24) 



This result suggests that the traditional decomposi- 
tion of the electron-phonon coupling parameter A into 
N(Ep)V ep is possibly misleading in the vicinity of band 
edges and van Hove singularities. 

In metallic systems without covalent bonding, there is 
a significant convergence problem in our method for cal- 
culating rj(E). For example in a simple metal such as 
aluminum, as shown in figure we cannot converge an 
answer with respect to either k-point sampling or super- 
cell size. The physical phenomenon causing this problem 
is that charge density changes in a metal such as the 
displacement of an atom lead to long range Friedel oscil- 
lations which are difficult to sample with the necessary 
precision. The magnitude of the electron-phonon cou- 
pling associated with this phenomenon is much smaller 
than the coupling due to strong covalent bonds. Any 
sampling errors associated with Friedel oscillations are 
negligible in systems with larger electron-phonon cou- 
pling contributions coming from covalent bonds. 

There is a large variation of electron-phonon coupling 
in the examples studied, a factor of five increase in the 
maximum of rj(E) in going from aluminum to silicon 
and another factor of five increase in carbon over silicon. 
These variations of magnitude can be approximately re- 
produced based on scaling arguments. The first approxi- 
mation we make is to assume that dEn^/dHi in Eq. 11' I'D 
is independent of n and k variations. This allows us to 
separate out the nk sum and equate it to the density 



of states per spin times the variation in the electronic 
energy, 



rk{E)f»N r (E) 



dE, 



'el 



(25) 



The sensitivity of electronic energy to atomic displace- 
ments is approximately related to the hydrostatic defor- 
mation potential - the sensitivity of electronic energy to 
change in the volume of the unit cell times the volume - 
divided by a length scale, in this case the bond length. In 
a free electron gas, the standard result for the hydro- 
static deformation potential is \Ep. Our simple model 
for rj(E) is thus 



■J \ 1 '■bona 



(26) 



The iV-p (E) is taken from figures |21 and [3] at the max- 
imum of rj(E) and from figure at Ep. The value of 
Ep in this model is taken to be the bandwidth of occu- 
pied valence states again taken from the figures: 22 eV 
for carbon, 15 eV for silicon, and 12 eV for aluminum. 
With the experimental bond lengths, we get rj max values 
of 27 eV/A 2 for carbon, 11 eV/A 2 for silicon, and 1.6 
eV/A 2 for aluminum. While these values are not quanti- 
tatively correct, they reasonable reproduce the trend of 
increasing electron-phonon coupling in these three mate- 
rials. The large variations in rj are just magnifications of 
the variations in bond lengths and strengths. 

The possibility of a high T c becomes more probable 
with increasing rji(E) and decreasing Mj for the atoms 
constituting a material. In going from silicon to carbon 
in a diamond structure we see that both these param- 
eters favor the lighter element. The lightest elements 
that can form strong covalent bonds under ambient con- 
ditions are hydrogen, beryllium, boron, carbon, nitrogen, 
oxygen, and fluorine. These are the elements that should 
make up the covalent framework central to the idea of 
a covalent metal superconductor. A hypothetical cova- 
lent metal superconductor need not exclude other ele- 
ments; they might play a role as dopants or be themselves 
bonded to one of the strong covalent elements. This en- 
compasses a broad class of materials, but two important 
considerations must guide and constrain the search for 
strong electron-phonon coupling. The features of rj(E) 
associated with covalent bonds are peaked in a relatively 
narrow range of energies. Bonds between different pairs 
of elements of different lengths and strengths will con- 
tribute to rj(E) at different energies. Only electrons at 
the Fermi energy can contribute to the superconducting 
state, therefore it is important to choose a material with 
bonds that primarily contribute their electron-phonon 
coupling to one energy. Additionally, if a material has 
bonds of different strengths, the weakest bonds will typ- 
ically have their bonding and anti-bonding states closest 
to the undoped Fermi level. If the Fermi level change 
required to bring rj(Ep) to the strong bonding (anti- 
bonding) states crosses weaker bonding (anti-bonding) 
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states, it might sufficiently deplete (fill) the weak bond 
enough to break it and destabilize the material. In some 
cases this problem can be averted by shifting the relative 
energy levels of bonds through charge transfer such as 
that which occurs with the a and 7r states of MgB2 26 . 

For a given material we can now roughly bound the 
maximum possible T c by rigidly shifting the Fermi level 
and using rj(E), Eq. and Eq. lfTU)l to tune to the 

largest T c . There are several caveats to this approach, 
the primary one being that synthesizing heavily doped 
materials can be very difficult. This difficulty arises from 
the presence of more energetically favorable structures or 
structural instability in the candidate material. Even if it 
is possible to physically modify a system and effectively 
shift its Fermi level, the modification might also change 
rj(E) from its initial value, for better or worse. An addi- 
tional problem arises in molecular crystals where weaken- 
ing intermolecular hopping and narrow bands can cause 
a divergence in both N(E) and rj(E) and a breakdown of 
Eliashberg theory. In the case of doped Ceo solids, it is 
believed^ 7 that non-adiabatic effects and a Mott insulator 
transition are required to fully account for some observed 
behavior. These are electronic effects beyond the scope 
of this paper that potentially bound T c by electronic en- 
ergy scales. More relevant to non-molecular crystals with 
a wider electronic bandwidth, we next examine a funda- 
mental instability of the phonon- mediated superconduct- 
ing state that bounds T c by a smaller, phonon energy 
scale. 



where i, j are the atom indices and a, (3 are the Cartesian 
coordinate indices. The first two terms are the classical 
electrostatic contribution between ions and a fixed elec- 
tronic charge background, while the third contains con- 
tributions coming from the electronic response. All Fermi 
surface contributions are contained within the third term 
of the dynamical matrix, and we rewrite it in a more sym- 
metric manner using the chain rule, 



dr 



dp e i(v) dV lon (r) 

9Rf OR? 
dV ion (r') dp e i(r) dV ion (r) 
dRf dV(r') dR p 



drdv'. 



(28) 



The middle term, dp/dV, is the electron susceptibility, 
usually denoted as x^ that quantifies the response of the 
electron density to a change in the background potential. 
If the contribution to \ due to the Fermi surface is small, 
then it can be related to the Fermi surface contribution 
to the irreducible susceptibility xo through perturbation 
theory^, 



FS 



(29) 



where e _1 is the inverse dielectric function of the mate- 
rial and is its Hermitian conjugate. The irreducible 
susceptibility can be expressed analytically in terms of 
electronic states, 



Xo(r,r') = 2 



nk.n'k' 



(7ik|r)(r|n'k')(n'k'|r')(r>k) 



IV. PHONON SOFTENING 

The phonon softening due to electron-phonon inter- 
actions can lead to a structural instability that sets an 
upper limit to the possible T c of a material. Phonon soft- 
ening is defined as a change in phonon frequencies from 
some well-defined bare frequency to the actual frequen- 
cies present in a material. The early research on phonon 
softening took the bare phonons to be of purely electro- 
static origin, as either coming from a screened ion-ion 
interaction^ or the ions and frozen valence electrons 7 -. A 
recent studyiS has included quantum effects in the defi- 
nition of the bare phonons, using metallic screening with 
an insulating electronic susceptibility. Since supercon- 
ductivity is the result of interactions between phonons 
and electrons near the Fermi surface, we will define the 
bare phonon frequencies here by isolating and removing 
all electronic contributions to phonon frequencies coming 
from the Fermi surface. 

To calculate the phonon spectrum of a material we 
must first consider its dynamical matrix, 



D 



a(3 



d 2 E 



d 2 E, t 



ion— ton 



dRfdR 



dRfdR? 



Pei(r) 



d 2 V lon (r) 
dRfdR? 



dr 



dp e i{r) dVio n (r) 
dRf dBS 



dr, (27) 



E n > 



E, 



nk 



(30) 



for Fermi occupations /„k- The susceptibility contains 
contributions from all pairs of electronic states, and to 
restrict these pairings to the vicinity of the Fermi sur- 
face we introduce energy cutoffs, Ep — AE < E < 
Ep + AE, to both sums over states. For sufficiently 
small AE, the second term in Eq. I|3U|I can be approx- 
imated by dfnk/dE, which can in turn be approximated 
by —5(Ep — E n k) for sufficiently low temperatures. With 
this as our expression for Ax$ s , we can now approximate 
the contribution of the Fermi surface electrons to the dy- 
namical matrix using Eq. (|29(l and Eq. I|28[) . 



AD 



af3,FS 



2 £ (n , k |^|g£| n '. k ') 



OR? 



x(n',k'\^^\n,k)S(E F -E ni 



OR 



:6 ((AE) 2 — (E F — £ n 'k') 5 



(31) 



Here, the screened self-consistent potentials are the result 
of the inverse dielectric function acting on the bare ionic 
potentials, 



dVsCF _ ^idVi on 



dRf 



dRf 



(32) 
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and 9 is the unit step function. 

As AE — * 0, the AD FS term vanishes. The Fermi 
surface does not have an infinitesimal effect on phonon 
frequencies, but rather it affects an infinitesimal number 
of phonons. For a phonon of wavevector q, we would 
calculate the softening due to the Fermi surface by tak- 
ing the expectation value of AD FS with respect to the 
phonon displacement vector. This leads to a substitution 
in Eq. I|3f(l of k' — > k+q due to momentum conservation. 
For sufficiently small q or in special cases of Fermi surface 
nesting, the condition \Ep — £Vi(k+ q )| < AE holds for k 
over a finite area of the Fermi surface and lead to finite 
phonon softening. As AE is made larger, more phonons 
will be softened by AD FS but the approximations made 
in deriving Eq. 1|31|1 eventually break down and the cor- 
respondence with electrons at the Fermi surface is lost. 
All the electrons in a material can play a role in softening 
the phonon frequencies from the ionic plasma frequencies, 
and the Fermi surface is not independent of the rest of the 
electronic structure. These additional phonon softening 
effects are not directly related to the Fermi surface and 
must be handled separately for different classes of ma- 
terials. For example, in metallic Lii-^BC-S^S there is 
an essentially constant softening for phonons from q = 
to |q| = 2kp. We imagine that for a hypothetical well- 
engineered superconductor, all material-specific phonon 
softening is designed so that it doesn't lead to any struc- 
tural instabilities, leaving only AD FS as the fundamen- 
tally limiting cause of phonon instability. 

With several more manipulations, the phonon soften- 
ing due to electron-phonon coupling can be directly re- 
lated back to a bound on T c . We first relate AD FS to 
the softening of an individual phonon mode with dis- 
placement vector £?, 

(Alu 2 ) fs = Y -S=AAf ' FS ^=- (33) 

Summing over all phonon modes at the gamma point, we 
have a complete basis of gamma point atomic displace- 
ments and can relate the phonon softening back to the 
total electron phonon coupling, 

£(Aa,?)™«2A<£5 2 ), (34) 

i 

assuming rji has been calculated using Eq. (|22[1 within 
the supercell approximation. This leads to a bound on 
T c by invoking the weak upper bound of Eq. to get 



T C <0T3 /£(Ao,?)™. (35) 

This phonon softening T c bound gets large as more 
phonons are softened and as the magnitude of the soft- 
ening increases. However, the total phonon softening is 
related to the total electron-phonon coupling strength, 
which is itself independent of the phonon spectrum. 



Phonon softening becomes the relevant bound on T c only 
when there is enough electron-phonon coupling in a ma- 
terial to drive a phonon mode to instability before the 
electron-phonon coupling strength is maximized. The 
identification of this one "Achilles' heel" mode and its 
subsequent destabilization are required to correctly in- 
terpret Eq. (|35|) . The degree to which a phonon mode 
can even be softened by the Fermi surface has to do with 
its overlap in Eq. (|33() with AD FS ; one has to evaluate 
these overlaps to determine the relative softening of vari- 
ous phonon modes. This "Achilles' heel" might not even 
be determined by softening due to electrons at the Fermi 
surface - other valence electrons and physical modifica- 
tions required to dope a material might both play an im- 
portant material-specific role in the softening of phonons. 
A caveat to this phonon softening argument is that it has 
been made based on the harmonic phonon approxima- 
tion which will cease to hold if some phonon modes are 
sufficiently softened. Anharmonic corrections to phonon 
frequencies can stabilize locally unstable modes, further 
delaying the onset of a true instability and raising the 
associated T c bound. 

A simple application of Eq. I|35|l is the case of di- 
amond. In diamond, the gamma point has three zero 
frequency acoustic modes and three equal frequency op- 
tical modes at the gamma point. The maximum phonon 
softening allowed in a stable structure is full softening 
to zero of the optical modes. The gamma point opti- 
cal modes have an experimentally measured frequency of 
165 meV— , which from Eq. i ffifll gives a T c bound of 430 
K. According to figure [3 the maximum T c possible with 
rigid hole doping is 290 K and with electron doping it is 
420 K. Within a rigid doping scheme, the theoretical limit 
of superconductivity in diamond is set by the amount of 
available electron-phonon coupling and not by a funda- 
mental lattice instability, due to the extreme stiffness of 
the undoped structure. 



V. T c BOUNDS OF COVALENT METALS 

Originating from the discovery of superconductivity in 
MgB2 in 2001—, there has been a growing interest in co- 
valent metals. One goal of this area of research is to 
determine, if possible, how to increase T c in already su- 
perconducting materials and how to drive insulating co- 
valent materials into a metallic, superconducting regime. 
The most prominent materials in these studies are two 
layered materials, metallic MgP>2 and semiconducting (at 
x = 0) Lii-^BC, and one isotropic material, boron-doped 
diamond B^Ci-^. The rji (-E)-based analysis derived in 
section IIIII provides a clear picture of the T c bounds in- 
herent in these materials. 

In MgB2 , neither hole-doping nor electron-doping have 
been successful in increasing T^. This has been ex- 
plained theoretically in terms of the magnitude of the 
phonon softening of the Ei g phonon mode as a func- 
tion of doping level. Because it is near to the strongly 
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FIG. 8: Density of states (dashed line) and rj(E) (solid line) 
for MgB2 . This calculation is performed on one unit cell with 
a 24 x 24 x 24 k-grid and 0.03 Rydberg gaussian smearing. 
Only in-plane motions of the boron atoms are considered; the 
out of plane boron motions and magnesium atoms do not 
contribute significantly to rj(E) in the valence bands. 



coupled regime, Awl, the variation of T c with doping 
should follow ?j(E), shown in figure |H1 This figure clearly 
shows that MgB2 is an optimally doped superconductor 
and that T c should dramatically decrease with electron- 
doping and gradually decrease with hole-doping, which 
accurately reproduces the results of Zhang et alifi. As 
for the magnitude of T c at optimal doping, our calcula- 
tion gives T c = 29 ~ 39 K using the Einstein spectrum T c 
formula of Eq. JTDJ} with fj,* = 0.1 - 0.15 and A = 0.75 
determined by setting \J (ui 2 ) — 76 meV, the anharmonic 
frequency of the E2 9 phonon modeS. This range of T c 
values encompasses the experimental value of 39 K for a 
physical range of /i* values. 

In LiBC, it has been predicted that superconductivity 
will occur upon uniform deintercalation of the lithium 
to the level of greater than ~ 10% and T c will plateau 
at ~ 50% deintercalation with values as low as 65 K 3 ^ 
and as high as 100 K 33 . The plot of r}(E) for LiBC in 
figure shows a magnitude of electron-phonon coupling 
twice that of MgB2 in figure [5] but less than half that of 
graphite in figure El Upon hole doping to the peak of 
rj(E) in the valence states, we predict a transition tem- 
perature in the range T c = 90 ~ 100 K using the Einstein 
spectrum T c formula of Eq. JTDJ) with ^* = 0.1 ~ 0.15 
and A = 1.6 determined by setting \J (u> 2 ) = 85 meV, 
the calculated frequency of the F,2 g phonon mode after 
doping^ (softened from 150 meV). The broad plateau 
in rj(E) for a wide range of doping levels corresponds to 
the T c plateau found in prior predictions 3 ^ 3 ^. An ex- 
perimental method for removing lithium or synthesizing 
lithium-deficient LiBC without degrading the structure 
has not yet been discovered and remains the barrier to 
fulfilling these predictions. 

In diamond, it has been observed that heavy boron 



FIG. 9: Density of states (dashed line) and rj(E) (solid line) 
for LiBC. DOS is scaled to half the actual (LiBC)2 unit cell to 
be comparable to the unit cells of AA-stacked graphite and 
MgB2. This calculation is performed on one unit cell with 
a 24 x 24 x 12 k-grid and 0.03 Rydberg gaussian smearing. 
Only in-plane motions of the boron and carbon atoms are 
considered; their out of plane motions and lithium atoms do 
not contribute significantly to fj(E) in the valence bands. In 
a single unit cell approximation, rji(E) is identical for boron 
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FIG. 10: Density of states (dashed line) and rj(E) (solid line) 
for zinc-blende boron carbide. This calculation is performed 
on one unit cell with a 24 x 24 x 24 k-grid and 0.03 Rydberg 
gaussian smearing. The equilibrium lattice constant is calcu- 
lated to be 3.72 A. In a single unit cell approximation, rji(E) 
is identical for boron and carbon. 



doping, up to the solubility limit, leads to superconduc- 
tivity with T c m 4 K in bulk diamond 34 and T c m 9 K in 
diamond thin films^. This doping level corresponds to 
~ 3% boron substitution and is still far from the optimal 
value of 50% boron substitution estimated from integrat- 
ing the DOS from Ep to the valence band peak of rj(E) 
in figure A small unit cell material corresponding to 
this optimal doping level is BC in a zinc-blende structure, 
which has been examined previous in the context of ma- 
terial strength 3 ^. We calculate rj{E) for zinc-blende BC 
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in figure 1101 and verify the optimal doping, but there is 
a substantial drop in the overall scale of electron-phonon 
coupling from the strong C-C bonds in diamond to the 
weaker B-C bonds. The scale of the electron-phonon cou- 
pling is comparable to the graphitic sheets of boron and 
carbon in LiBC, but the peak in the zinc-blende struc- 
ture is narrower and ~ 20% higher. This material is 
stable, unlike BC in a graphitic structure, with a dra- 
matic softening of the gamma point optical mode from 
165 meV in diamond to a calculated value of 55 meV 
in BC. This phonon calculation is performed using the 
frozen phonon method with anharmonic corrections; the 
total energy as a function of phonon displacement is fit 
to a 6 th order polynomial and the frequency is calcu- 
lated as the energy difference between the first excited 
state and ground state of the associated anharmonic os- 
cillator. From the calculated phonon frequency, we get 
A = 3.6 and using the Einstein spectrum T c formula of 
Eq. JTHJ with ijl* = 0.1 ~ 0.15, we predict the tran- 
sition temperature of this material to be in the range 
T c = 140 ~ 160 K. The barrier to reaching T c values of 
this magnitude in diamond is the difficulty of including 
boron into the material beyond the solubility limit and 
the formation of boron-based impurities that don't act 
as acceptors^!. 

A recent study^ has examined the possibility of 
electron-doping in diamond, which they conclude might 
lead to larger T c values than hole-doping. This study was 
based on small doping levels, less than 6%, and we con- 
tinue this argument to the doping limits. The rj(E) peak 
in the conduction band of diamond is more than twice 
as large as the valence peak, but requires significantly 
more doping to reach. Electron-doping is presumed to be 
achieved by nitrogen substitution and a material repre- 
sentative of 50% doping is CN in a zinc-blende structure, 
which is known to be structurally unstable^. The rj{Ep) 
of zinc-blende CN is smaller than in zinc-blende BC and 
still far from the peak in the conduction band. We con- 
clude from this that hole doping is the best prospect for 
exploring the limits of T c in superconducting diamond. 

VI. CONCLUSIONS 

The phonon frequency being a bound on T c is an ob- 
vious conclusion from the weak coupling limit of con- 
ventional electron-phonon superconductivity. It takes a 
great deal more effort to show that this is still a bound in 
the strong coupling limit where T c at first glance appears 



not to have a bound. This bound is caused by the soften- 
ing of phonon modes to instability due to strong electron- 
phonon coupling. In a simplified Eliashberg framework, 
the variation with doping and hypothetical maximum of 
electron-phonon coupling can be studied through a sin- 
gle efficiently calculable metric, rj(E). In a material with 
strong covalent bonds, such as diamond, we can see that 
the T c bound due to softening cannot be reached because 
there is not enough electron-phonon coupling strength 
available to soften the phonons to instability. In such a 
material, an upper bound on T c is set by the maximum 
available electron-phonon coupling. 

Of the materials studied, only MgB2 has maximized 
its electron-phonon coupling within an as-yet experimen- 
tally realizable material. Careful materials design and 
synthesis will be required to increase the maximum avail- 
able electron-phonon coupling and approach the T c limits 
set by structural stability. For some materials contain- 
ing carbon-carbon and equally strong bonds and light 
atomic masses, this T c limit is higher than room temper- 
ature. These hypothetical materials are likely to need 
weak links to enhance N(Ep) and therefore rj(Ep) and 
interstitial regions to contain donor or acceptor atoms or 
molecules that don't play a role in the stability of the 
lattice and can be removed or added to chemically dope 
the covalent network into a metallic regime. Molecular 
crystals fit this description, but intermolecular hopping is 
often too weak and prevents conduction. Polymerization 
of some of the constituent molecules might serve to in- 
crease hopping and enhance conduction. These criteria 
constitute a broad class of hypothetical materials, and 
further theoretical progress can be directed and stimu- 
lated by the experimental realization of such materials - 
much as the recent interest in phonon-mediated super- 
conductivity in covalent metals has been motivated by 
the discovery of superconductivity in MgB2. 
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